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Abstract 

A reliable rotor aeroelastlc analysis opera- 
tional at Hughes Helicopters, Inc. that correctly 
predicts the vibration levels for a helicopter Is 
utilized for the present study to test various 
unsteady aerodynamics models with the objective 
of Improving the correlation between test and 
theory. This analysis called Rotor Aeroelastlc 
Vibration (RAVIB) computer program Is based on a 
frequency domain forced response analysis which 
utilizes the transfer matrix techniques to model 
hell cop ter/rotor dynamic systems of varying 
degrees of complexity. The analysis Is a n on- 
modal analysis and It Includes effects of periodic 
coefficients for the forward flight conditions. 

The first new aerodynamics model Incorporated In 
the analysis was based on the current state-of-art 
of unsteady aerodynamics. The results based on 
this aerodynamics model for the AH-1G helicopter 
rotor were compared with the flight test data dur- 
ing high speed operation and they Indicated a 
reasonably good correlation for the beamwlse and 
chordwlse blade bending moments, but for torsional 
moments the correlation was poor. As a result, a 
new aerodynamics model based on unstalled synthe- 
sized data derived from the large amplitude oscil- 
lating airfoil experiments was developed and 
tested with RAVIB analysis. The results Indicate 
a significant Improvement In the correlation for 
the torsional moments. 

Introduction 

The rotor aeroelastlc stability and response 
analyses (e.g.. Reference 1) Invariably Involve 
computation of airloads through use of unsteady 
strip theories. To a significant extent, accuracy 
of analytically predicted parameters, such as aero- 
dynamic damping, blade oscillatory bending and tor- 
sional moments and hub vibratory loads, depends on 
the correctness of the unsteady aerodynamics model 
utilized In the analysis. Recently there has been 
significant effort (References 2-5) to develop new 
unsteady aerodynamics models for the rotary-wing 
applications with a correct emphasis on Incorpora- 
tion of precise nature of wind and rotor blade 
motions. But It has not been fully demonstrated 
(through correlation with test data) the kind of 
Improvements we can expect by utilizing these 
more sophisticated unsteady aerodynamic models, 
iiie present study was undertaken to test various 
unsteady aerodynamics models with the primary 
objective of Improving the correlation between the 
flight test data and theory. 

Background 

Most of the unsteady aerodynamics models 
involve extension of Theoderson's flat plate 
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theory (Reference 6) to account for the airfoil 
pitch-plunge motion and also to represent the 
effects of time- varying free stream velocity. The 
airfoil shape and thickness effects are normally 
accounted through replacement of the flat plate 
lift curve slope (2ir) by an appropriate airfoil 
static lift curve slope. Blade section unsteady 
lift and pitching moment coefficients normally 
Involve a lift deficiency function C*(k); where k 
is the reduced frequency. For ro^ry wing applica- 
tions, where blade section experl *rces arbitrary 
motions, use of reduced frequency Is highly Inap- 
propriate. Therefore, In rotor aeroelastlc analy- 
ses, It Is practical to model unsteady aerodynamics 
of airfoil arbitrary motions preferably In both 
Laplace and time domains. This Is simply achieved 
by converting the generalized lift deficiency func- 
tion Into Pade' form (e.g.. Reference 7) where It 
1$ directly available In Laplace domain. For the 
time-domain applications, Pade' form can be easily 
converted Into highly practical Indlclal formula- 
tions, such as indlclal form of Wagner function for 
the oscillating flat plate. 

Despite availability of the above-mentioned 
methodology for unsteady strip theory, most of the 
rotary-wing dynamlclsts find It convenient to just 
replace the lift deficiency function by a constant 
number. Conventionally a value of 0.8 to 1.0 for 
C*(k) Is used. To date, an uncertainty exists 
regarding the benefits of these more sophisticated 
ursteady aerodynamics theories. However, the pre- 
sent study clearly indicates that conventional 
aerodynamics models are highly Inaccurate In pre- 
dicting certain blade section airloads, such as 
unsteady pitching moment. More specifically, a 
new but conventional unsteady aerodynamics model 
(Including aerodynamic spring-damper matrices) 
based on current state-of-art of the rotary-wing 
aerodynamics (Reference 8) was developed and It 
was Incorporated Into a reliable rotor aeroelastlc 
analysis operational at Hughes Helicopters, Inc. 
(Reference 1). The computed results (airloads and 
oscillatory blade bending and torsional moments) 
for AH-1G rotor blade were compared with the avail- 
able flight test data corresponding to high speed 
flight conditions. The correlations obtained for 
the beamwlse and chordwlse bending moments were 
good but for torsional moments, the correlation 
was poor. As a result, a new unsteady aero- 
dynamics model based on unstalled oscillating air- 
foil test data was developed and Incorporated In 
the analysis. The results Indicate a significant 
Improvement In the correlation of computed oscil- 
latory torsional moments with test dtta. 

Description of Aeroelastlc Analysis 

The aeroelastlc analysis utilized for the 
present study Is called RAVIB (Rotor Aeroelastlc 
Vibration) computer program and It Is a modified 
and Improved version of a computer program 
(References 1 and 9) originally developed by 
Rochester Applied Science Associates. Briefly, 
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the RAVIB computer program is based on a frequency 
domain forced response analysis which utilizes 
the transfer matrix techniques co model helicopter/ 
rotor dynamic systems of varying degrees of com- 
plexity. The analysis is a non-model analysis 
and it involves application of a standard matrix 
process in which the transfer matrices associated 
with successive characteristics of the modeled 
blade are combined to form the transfer or associ- 
ate matrix relating the shears, moments, slopes, 
and deflections at a position on the blade to 
those occurring at the tip of the blade. The 
blade mass and structure characteristics are 
represented in a lumped parameter form. The analy- 
sis Includes aerodynamic interharmonic blade 
coupling (periodic coefficients) and interhar- 
monic coupling due to fuselage motion. Only a 
few features of the analysis will be breifly dis- 
cussed here. More details are given in Refer- 
ence 1. Briefly, rotor/helicopter model consists 
of a rotor system with flexible blades which may 
be articulated, gimballed, teetering or hinge- 
less type. The rotor may be connected to fuselage 
through a fixed system rotor support consisting of 
gearbox with roll and pitch flexibility and a 
flexible drive shaft. Fuselage may be modeled as 
a flexible beam (similar to blade model but non- 
rotating) or fuselage effects may be represented 
by a hub impedance matrices. The program has 
capability to model a detailed swashplate-type 
control system. The basic rotor blade structure 
is represented by a lumped parameter model in 
which the blade is subdivided into a finite num- 
ber of blade sections. Each blade section may 
have (see Fig. 1) arbitrary orientation and chord- 
wise location of shear center, arbitrary spanwise 
distribution of mass and Inertias, twist, chord, 
mass C.G. location, bending and torsional stiff- 
nesses, chordwise aerodynamic center location. 

The aerodynamic effects Include aerodynamic 
inertia, damping and spring rates which vary azi- 
muthal ly (periodic coefficients) in forward 
flight. Radial and azimuthal variations of wake 
Induced velocities may be included Also, 
deformed free-wake effects on helicopter rotor 
system dynamic response may be included in the 
analysis In an iterative procedure which couples 
RASA free-wake analysis (Reference 10) to the 
blade motions. 
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Most significant of these modification, are listed 
below: 

1) An Iterative procedure was Incorporated 
in the analysis to obtain the compatible steady 
state elastic deflections of the blade. A number 
of sources of blade vibratory excitations (for 
example, blade aerodynamic pitching moment) vary 
in a highly nonlinear fashion with the steady 
state elastic deflection of the blade. Inclusion 
of this procedure successfully minimizes the error 
due to these nonlinearities. 

2) The RAVIB analysis was modified to couple 
the blade root motions with fuselage through a hub 
impedance matrix and thus to account for the fuse- 
lage motion effects. It is assumed the impedance 
matrix can be conveniently obtained by exercising 
NASTRAN. 

3) Unsteady aerodynamic effects of higher 
harmonic control (HHC) Inputs were incorporated 
in the analysis. 

4) A new unsteady aerodynamic model based 
on the current state-of-art aerodynamics was 
developed and it was successfully incorporated 
in the analysis. Further discussions of this 
unsteady aerodynamics modi 1 is provided next. 

Description of Unsteady Aerodynamics Model 

The aerodynamic forces acting on a blade 
section in time-domain, as it goes around azimuth, 
can be represented in matrix form by the follow- 
ing relationships: 


{F A } = {F Ao } + DV(V> ♦ [C A J {$} + [k A J {q} (1) 

where {FAo> is a six-component vector of aerody- 
namic forces and moments in an appropriate coordi- 
nate system due to all known motions of blade and 
wind. The matrices [M/\] , [Cfl] and [k/\] represent 
mass, damping and spring rates respectively due 
to aerodynamic forces. The vector {q} is a six- 
component state vector of unknown deflections and 
rotations. In general, for steady state heli- 
copter flight conditions, the matrices [M^], 

[Ca] and [k/^] vary periodically around the azimuth. 
Furthermore, if Theoderson-type unsteady aero- 
dynamics are utilized, these matrices contain 
lift deficiency function C'(k) in one form ,r 
another. Because it is not possible to precisely 
describe the reduced frequency k for the blade 
section, some approximations are required to evalu- 
ate C * (k) during the computations. The most con 
mon approximation involves assuming a constant 
value of 0.75 to 1.0 for C'(k), The other common 
procedure Involves transforming C'(k) into 
indiclal form (Reference 7) and computing the air- 
loads in time-domain. This procedure is highly 
practical for computing [Fa 0 ], where time-history 
of aerodynamic angle of attack is completely 
known. The most commonly used indicia! form for 
C(k) Is the Wagner functions (Reference 11) 
derived from the flat plate theory. 


Recent Improveme r 

As mentioned earlier, the RAVIB analysis has 
been developed by modifying an existing analysis 
developed by Rochester Applied Science Associates 
(Reference 1). A number of these modifications 
or Improvements were absolutely necessary for 
obtaining a good correlation between predicted 
results and test data. A few other modifications 
were carried out to enhance its capabilities. 
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Frequency Domain Formulation for Airloads 

The present study requires development of 
aerodynamic transfer matrices in Laplace form 
thic can be used in the frequency-domain analysis 
of Reference 1. Except fur a small magnitude 
terms Involving the function [1-C'(k)]» each of 
these periodically varying matrices [M A ], [Ca] and 
[ka] can be expanded In a Fourier series form. 

For example, 




[C A ] * I [C n 3e innt (2 

n=-°° 

Similarly, the periodic motion of the blade 
around the azimuth can be represented as follows: 


W - I l\l (3) 

k=-® 


Thus the aerodynamic forces generated by elastic 
deflections (F^q) involve multiplication of two 
infinite series resulting in interharmonic 
coupling as follows: 



*( jjg ■“)(* jjv •**) 


(4) 


Thus, due to periodic coefficients, for example 
during high speed forward flight, a significant 
amount of interharmonic coupling occurs. 


The Laplace transform is used to transform 
the differential equations into a set of algebraic 
equations as follows: 


!U ■ 2 [< -ikft-inn) 2 [M n ] + (-ikn - inti) [C p ] 
n=-“ 


+ |q k+n | (5) 

for k * -«a, It 2, 3, 4 ...» 


Here vector {q^} represents Laplace transform of 
k*" component of the deflection vector {q}. In 
principle, Equations 5 can be used to solve any 
number of harmonics simultaneously. But In prac- 
tice, It Is sufficient to truncate the summation 
(n 9 -1, 0, +1) for each value of k. For example, 
if aerodynamic transfer matrix corresponding to 
4 per rev (k-4) response Is desired, the blade 
harmonic motions at 3, 4 and 5 per rev 
(qo Qfj) have to be simultaneously computed 
under this procedure. The error Involved is pre- 
sumed to be small due to the exclusion of higher 
order Interharmonic coupling (n>2). Moreover, If 
desired, once all the desired harmonics (Np) 
have been computed (q^, k ■ 1, 2,..,N p ) by above 


procedure, more accurate airloads can be computed 
during a final "pass" using Equations 5 for full 
N p values of n. Thus, an iterative procedure will 
reduce the error due to the truncation to a mini- 
mum. There are a number of additional factors 
which may make it necessary to follow the iterative 
procedure in order to obtain more accurate results. 
These are discussed next. 

Possible Reasons for Iterative Procedure 

1) As mentioned in the previous section, the 
periodic coefficient matrices involve terms pro- 
portional to [l-C'(k)], which cannot be conve- 
niently accounted for under the present procedure; 
unless, of course, C’(k) is given a consent value. 
Even though these terms are small, they can be 
appropriately included if an iterative procedure 

is followed. During each iteration, vector { q) 
corresponding to previous iteration can be used 
and [l-C'(k)] can be replaced by an indicial form 
equivalent (for example indicial form of Wagner 
function, Reference 5, Page 15). Thus, the result- 
ing time-domain forces corresponding to [l-C'(k)] 
can be computed and these can be included in the 
vector {Fao) of the Equations 1. 

2) If the elastic deflections of the blade 
are significantly large, such as large torsional 
deflections during dynamic stall, the matrices in 
Equation 1 may be in error due to presence of 
significant nonlinearities. But, if an iterative 
procedure is followed, the aerodynamic forces 
corresponding to large blade elastic deflections 
(estimated from previous iteration) can be directly 
included into forcing vector f F/\ 0 } of Equation 1. 
This procedure will minimize the errors due to the 
nonlinearities . 

3) For correct blade dynamic response, it is 
necessary to use nonuniform induced velocity dis- 
tribution over the rotor disk. This nonuniform 
induced velocity distribution is normally com- 
puted by using a free-wake analysis (Reference 10, 
fur present analysis), wherein strength of wake 
vortices depend on blade dynamic response. Thus, 
deformed free-wake effects on rotor dynamic 
response can be correctly accounted only by fol- 
lowing a solution method which Involves iterative 
procedure that couples the free-wake circulations 
to the blade motions. 


Application of New Conventional Unsteady Aero- 
dynamics Model 

The existing unsteady aerodynamics model in 
Reference 1 analysis was found to be highly inac- 
curate and limited in scope. For example, it 
assumed the aerodynamic mass matrix in 
Equation 1 to be equal to zero. A new unsteady 
aerodynamics model based on current state-of-art 
was developed and Incorporated in the RAVIB analy- 
sis (Reference 1). The basic equations utilized 
were similar to the ones given In Reference 6 
(Section 11-8 on Page 596), The effects of 
radial flow, lift deficiency function, dynamic 
inflow (optional) were appropriately Incorporated 
in the model. The development essentially 
involved accurate formulation and programming of 
the aerodynamic forcing function {Fao' and the 
matrices [MaL and [k^] (see Equation 1) 
in the RAVID analysis. The resulting analysis 
was utilized to carry out a correlation of com- 
puted blade aeroelastlc airloads, bending and 
torsional moments with flight test data with an 
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objective of establishing the airloads model. 

During an earlier study (Reference 12), the fre- 
quencies and mode shapes obtained from the Ref- 
erence 1 analysis were compared with those from 
other independent analyses and the two results 
compared very well. Thus, the analysis to be 
utilized for the present study is a validated 
computer program that is believed to represent 
dynamics of helicopter blade quite accurately. 

Application of Analysis 

The RAVIB analysis has been used to compute 
blade airloads, blade bending and torsional 
moments for the AH1-G helicopter flying in 
steady state flight conditions. For the AH1-G 
helicopter, flight test data are readily avail- 
able (Reference 13) through DATAMAP (Reference 
14). Therefore, it was found convenient to use 
the Reference 13 flight test data for the present 
correlation study. 

The AH1-G blade structural, geometric and 
aerodynamic characteristics were obtained from 
References 13 and 15. These blade characteristics 
may not correspond exactly to the AH1-G helicopter 
blade used in the flight tests. Because the main 
objective of the present study was to demonstrate 
relative improvements in the correlation between 
test and theory through the use of better unsteady 
aerodynamics model, an approximate but representa- 
tive model of the AH-1G blade was considered to be 
adequate. The blade was modeled by approximately 
twenty nonuniform elastic segments. The effects 
of parameters such as control system stiffness, 
undersling and drive shaft torsional flexibility 
were appropriately included in the model. The 
effects of hub impedance and drive system damping 
were neglected in the present computations. 

Results with Conventional Unsteady Aerodynamics 

The conventional strip theory was used for 
computation of airloads. The blade was represented 
by nine nonuniform strips. At the center of each 
strip (aerodynamic load application point), the 
nonuniform azimuthal distribution of induced 
velocities was computed by utilizing a rotor wake 
analysis (modified version of Reference 10). The 
computations for the AH-1G helicopter rotor were 
carried out at two level flight conditions (114 
knots and 142 knots) for which the flight test 
data were available (flight numbers 610 and 
614 in Reference 13). The gross weight for both 
the flight conditions was 8,300 pounds with 
Cj/o ® 0.006. 

A correlation between computed and measured 
performance parameters is shown in Fig. 2. As 
the results indicate, the predicted values of 
shaft horsepower and collective are within five 
HP and within half a degree respectively of the 
test values (Fig. 2). For a radial station cor- 
responding to r/p * 0.75, correlation between 
calculated and test airloads is shown in Fig. 3a 
for a representative flight condition (114 knots). 
Even though there are significant differences 
between test and theory when the blade is In cer- 
tain segments (e.g., near ty ■ 270) of rotor disk, 
overall correlation between the two Is reasonably 
good. The similar results for the 142 knots 
fMgnt condition are shown in Fig. 3b. 
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Fig. 2 Comparison of measured and computed level 
flight performance parameters, AH-1G 
helicopter, 8300 pounds, GW 

Figs. 4a and 4b show the blade oscillatory 
beamwlse bending moment correlation between the 
test and theory based on the conventional aerody- 
namics model. For 114 knots flight condition 
(Fig. 4a) the theory predicts peak to peak value 
quite accurately, but for 142 knots case (Fig. 4b), 
the theory underpredicts peak to peak value by 
30 percent. Time histories of measured and com- 
puted chordwlse bending moments for the same 
flight conditions are compared In Figs. 5a and 
5b. From the results shown In these figures, It 
Is seen that very good agreement between test and 
theory has been obtained over complete range of 
azimuth. 

The correlations of the test blade torsional 
moments with the predicted results are shown In 
Figs. 6a and 6b. The peak to peak variation In 
torsional moment Is highly underp^edlcted by the 
theory based on conventional aerodynamics model. 

A thorough analysis of these results Indicated that 
a new aerodynamics model was needed to Improve the 
correlation between computed and test torsional 
moments . 

Description of New Unsteady Aerodynamics Model 

An analysis of results from two-dimensional 
experiments (e.g., Reference 16) Involving large 
amplitude oscillations of airfoils (under unstalled 
conditions) Indicates that conventional aero- 
dynamics (based on Theoderson's theory) Is unable 
to predict the unsteady aerodynamic characteristics 
of the airfoil. This Is partly due to the fact 
that Theoderson's flat plate theory is based on 
small amplitude oscillations. During high speed 
forward flight, helicopter blade sections are 
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Fig. 5 Comparison between measured and computed 
oscillatory chordwlse bending moment, 
r/R - 0.39, AH-1G blade 

Fig. 6 Comparison between measured and computed 
oscillatory torsion moment, 
r/R - 0.30, AH-1G blade 
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expected to go through significantly large ampli- 
tude oscillatory changes In angle of attack. 

Thus, It was found necessary to Include the 
effects of finite amplitude oscillations In the 
new unsteady aerodynamics model. More specifi- 
cally, the new unsteady aerodynamics model util- 
izes empirical o\- synthesized data. The syn- 
thesized data are derived from the test data; 
by curve-fitting the appropriate analytical 
expressions to the measured unsteady airfoil 
characteristics, obtained from the oscillating 
airfoil experiments. The present method used 
the analytical expressions established in Ref- 
erence 5. For example, as described in Ref- 
erence 5, the unsteady lift coefficient, under 
unstalled conditions, is represented by 

C Lu = C Ls (a) + Q 1 A + Q 2 a w + Q 3 a + Q 4“ 2 < 6) 


Here Cl s is the static lift coefficient; (A, a w ,a) 
are the instantaneous values of dynamic parameters 
(see Appendix for details) and Qj, Q ? , Q3 and 
Q4 are the empirical parameters or the synthesized 
data. Because the empirical parameters Qj through 
Q4 are based on real airfoils executing large amp- 
litude oscillations, the Equation 6 correctly 
simulates general unsteady lift characteristics 
of a helicopter blade section. Furthermore, It 
should be noted that lift deficiency function 
effects are represented in Equation 6 through the 
decay parameter a w , which Is derived from a modi- 
fied Wagner function (Reference 5). In fact, 
in future, th^ parameters such ar Qj through Q 4 
can be obtained from some reliable analyses. 

Finally, in unsteady aerodynamics the ques- 
tion of what constitutes a large amplitude 
motion, depends to some extent on the magnitude 
of the Mach number. An amplitude of one degree 
in transonic flow is considered high amplitude; 
whereas at low subsonic Mach numbers, a three 
degree oscillation may be considered a small amp- 
litude motion. Studies are continuing to estab- 
lish these criteria. 

Next, the results based on this synthesized 
unsteady aerodynamics model are discussed. 

Results with Synthesized Unsteady Aerodynamics 

Equation 6 describes the unsteady lift 
coefficient, Cl u , of an airfoil in the time 
domain. Similar equation for unsteady pitching 
moment coefficient Is given In Reference 5. 

These two equations were Incorporated in the 
RAVIB analysis. The modified analysis was 
utilized to recompute the earlier results 
obtained from the conventional unsteady aero- 
dynamics theory. Thus, the various improvements 
in the correlation between flight test data and 
the analytical computations can be systematically 
demon stra ted. 

First, the variations in the predicted air- 
loads with aerodynamics model are shown In Figs, 

7 through 10. Fig. 7 has beer repeated from 
Fig. 3, but with the addition of computed normal 
force coefficient based on synthesized data. 

The differences between the two analytical 
results are small and these differences are 
mainly confined to the retreating blade region 
of the rotor disk. To further Illustrate the 


closeness of the computed Ci from the two aero- 
dynamic models, Figs. 8a and 8b show similar 
results obtained at one more radial station cor- 
responding to r/ R * 0.86. A close analysis of the 
results, however, indicates a small but signifi- 
cant improvement due to the use of synthesized 
data, but only for 142 knots flight condition 
(Figs. 7b and 8b). Thus, improvements in the 
correlation of beamwise bending moments, if any, 
can be expected only for 142 knots flight condi- 
tion. For completeness, Figs. 9 and 10 show the 
corresponding correlation of chordwise force 
coefficients for the two radial stations (0.75 
and 0.86). Because the results for the chordwise 
force coefficient from the two aerodynamic models 
differ only slightly, no significant variations 
in correlation for the chordwise bending moments 
is expected with change in the aerodynamics model. 
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Fig, 7 Correlation between calculated and test 
airloads, AH-1G blade, r/R = 0.75 
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Fig. 10 Correlation between calculated and 
test airloads, AH-1G blade, r/R = 0.864 

Bending Moment Correlations 

Figs. 11 and 12 show the comparison between 
the test and the predicted beanwise bending 
moment values described for two different radial 
stations, r/ R * 0.39 and r/ R = 0.80 respectively. 
An analysis of these figures indicates the two 
aerodynamics models predict similar results at 
114 knots flight condition (Figs. 11a and 12a). 
However, at 142 knots flight conditions, the 
predicted results based on synthesized data 
(Figs, lib and 12b) seem to compare better with 
the test data than those based on conventional 
aerodynamics. 

The similar variation of oscillatory chord- 
wise bending moments for various aerodynamics 
models are shown In Figs. 13 and 14, correspond- 
ing to two radial stations r/ R ■ 0.39 and 
r/ R * 0.8 respectively. Unlike beamwlse bending 
moments, the computed time histories of the chord- 
wise bending moments do not show significant vari- 
ation with the unsteady aerodynamic model. It 
should be remembered, however, that the main pur- 
pose for developing the present synthesized data 
model was to Improve the correlation between pre- 
dicted torsional moments and the test data. 
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Fig. 11 Comparison between measured and computed 
oscillatory beamwlse bending moment, 
r/ R = 0.39, AH-1G blade 


Torsional Moment Correlation 


Measured and computed torsional moments for 
114 knots flight condition are compared In Figs. 
15a and 16a for two different radial stations, 
r/R - 0.3 and r/R 3 0.5 respectively. From these 
figures, It is seen that very good agreement 
between test and new theory based on synthesized 
data has been obtained over the complete range 
of azimuth. The computed torsional moments based 
on conventional aerodynamics do not correlate well 
with the test data. The similar correlations are 
obtained at 142 flight conditions as Indicated by 
the results shown In Figs. 15b and 16b. Thus, the 
results shown in Figs. 15 and 16 Imply that for 
correct computations of pitch link loads, the 
unsteady aerodynamics model has to oe based on 
large amplitude Incidence oscillations Inherently 
encountered by a helicopter blade section. 


Finally, It should be mentioned that use of 
synthesized data will make it easy to extend the 
airloads computations In the stall region. As 
established In Reference 5, the unsteady aero- 
dynamics characteristics of blade section during 
dynamic stall are easily obtained by adding more 
terms to the right-hand side of Equation 6. 
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Concluding Remarks 

Rotor Aeroelast' *' ration (RAVI6) com- 
puter program is u islve, elegant and 

efficient analysis wnsch seems to predict blade 
oscillatory loads rccuy; ably well. The conven- 
tional unsteady aerotiyna;ni ;s mod^l based on small 
amplitude airfoil oscl rations seem to highly 
underpredict helicopter H iade pitching moments. 


Use of a new a*?r:j ics model based on unstalled 
synthesized data, "’e :>~'d from the large amplitude 
oscillating airfoils «. ts, signlflcalntly Improves 
the correlation of + ; «.. imputed blade torsional 
moments with the .ata. Furthermore, this 
ne ’ aerodynamics is such that It can be 
easily extended w ^te the critical pitch- 
link loads of * c i pirating under dynamic stall 
conditions. 



Fig. 12 Comparison between measured and 
compute! oscillatory beamwlse bending 
moment, r/ R - 0.80, AH- 16 blade 


Fig. 13 Comparison between measured and computed 
oscillatory chordwlse bending moment, 
r/ R ■ 0.3*, AH- 16 blade 
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Fig 16 Comparison between measured and 
computed oscillatory torsion 
moment, r/o * 0.50, 

AH-1G blade 

APPENDIX 


(b) 



The dynamic parameters utilized in the 
synthesize! data aerodynamics model are: 1) the 

instantaneous angle of attack, a\ 2) the nondi- 
menslonal >1 tch rate A; and 3) the decay para- 
meter a w , \hich accounts for the time history 
effects of the change in a, and Is based upon the 
Wagner function. 


Because the motion of a helicopter blade is 
not known a priori , the blade section dynamic 
parameters are evaluated numerically in a step- 
wise manne** by utilizing the following recursive 
relationships at step n. 


o - 6 _ + 4 > 
n n n 


(A-l) 


A n - k9) n /(*>n 


(A-2) 


+ (1.54» n - 2. ()♦„_! + 0.5* n . 2 )/tAs) n 


- X n + Y n 


(A-3) 


where 


-0.0455(1-M 2 ) (As) 

Vl e 


+ 0.165 (•„ - 

(A-4) 

-0.30(1-M 2 ) (As) n 

Vl e 


+ 0.335 (o n - 

(A-5) 

2U 

(as) = (a4>) 

(A-6) 


Here M is azimuthal steps ize, si is rotor 
speed, c is chord length, and U n is tangential 
velocity component. 

The instantaenous angle of attack, o n , is 
described in the tip-path-plane system, G n and <fr n 
being the pitch angle and inflow angle, respec- 
tively. It should be noted that the time deriva- 
tive of pitch angle in Eq. (A-2), (A6/A$), may be 
computed analytically from the known cyclic or 
harmonic Inputs, while the time derivative of <i> 
has to be computed by the backward difference 
scheme. The derivation of Eqs. (A-l thru A-6) for 
a w is described in Reference 5. 
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DISCUSSION 
Paper No. 8 

DEVELOPMENT OF AN UNSTEADY AERODYNAMICS MODEL TO IMPROVE CORRELATION 
OF COMPUTED BLADE STRESSES WITH TEST DATA 
Santu T. Cangwani 


Jing Yen, Bell Helicopter : I have two questions. Number one: Does your unsteady aerodynamics 

model include lift, drag, and moment — everything? 

Cangwani : The conventional model, you know, it includes only lift coefficient and pitching 

coefficient, and drag is computed mainly from the lift. In this model you still use the static 
data and most of the unsteady drag comes from the unsteady lift or from the static drag. 

Yen : I know from my experience that the conventional or well-known unsteady aerodynamics model 

can give you close to good correlation— the moments and everything else, but is very poor in 
power, and the drag. Your unsteady aerodynamics model here includes everything, including the 
drag? 

Cangwani : I showed you the correlation on horsepower. 

Yen : Yes, it was very good. 

Gangwani : Edgewise stress correlation is also very good so I don't see any problem in computing 

the edgewise pulses really. 

Yen : My next question is wiere did your synthesized data come from? 

Gangwani : The synthesized data I obtained as I told you from the oscillating airfoil test. Now 

your question is how did I get [the data] for a particular airfoil? 

Yen : That's right. 

Gangwani : Well, I sort of took it from other airfoils since your airfoil was a symmetric air- 

foil. So most of the data corresponds to a [NACA] 0012. It's not an absolute correlation, it's 
Just qualitative to demonstrate that we do need sor* different aerodynamics models to compute 
the aerodynamic pitching moments. 

Bob Blackwell , Sikorsky Aircraft : I want to ask one question. I wasn't completely clear about 

whether this model truly handles stall conditions or could it be extended so that it would? The 
second part of that is were the conditions shown at 142 knots and so forth, at that Cj/o, did 
th*»y represent conditions that really would design the control system? In other words were 
those loads the high loads for the system? 

Gangwani : The highest data available are at 142 knots and there was no stall. The stall 

results that I showed were Just computational results at Gf/o of 0.1. There was no test data 
at that point. They were Just computational data. Most of the test data are at . C^/o much 
lower than that. 
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